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We find that the polarization field, B x , obtained by scaling the weak-parallel-field magnetore- 
sistance at different electron densities in a dilute two-dimensional electron system in (111) silicon, 
corresponds to the spin susceptibility that grows strongly at low densities. The polarization field, 
-B sat , determined by resistance saturation, turns out to deviate to lower values than B x with in- 
creasing electron density, which can be explained by filling of the upper electron subbands in the 
fully spin-polarized regime. 
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Much interest has been attracted recently to the 
anomalous properties of strongly interacting two- 
dimensional (2D) electron systems. The interaction 
strength is characterized by the Wigner-Seitz radius, 
r s = l/(7rrt s ) 1 / 2 as (where n s is the electron density 
and ob is the effective Bohr radius in semiconductor), 
which is equal in the single- valley case to the ratio of the 
Coulomb and the Fermi energies. While at high electron 
densities (r s < 1) conventional Fermi- liquid behavior is 
established, at very low electron densities (r= 3> 1) for- 
mation of the Wigner crystal is expected [H, 0, El • It was 
not until recently that qualitative deviations from the 
weakly-interacting Fermi liquid behavior (in particular, 
the drastic increase of the effective electron mass with 
decreasing electron density) have been found in strongly 
correlated 2D electron systems (r s > 10) 0, [E@|- 

The strongest many-body effects have been observed in 
(lOO)-silicon metal-oxide-semiconductor field-effect tran- 
sistors (MOSFETs). This has particularly stimulated 
a new series of experiments in (lll)-silicon MOSFETs 
@, H, OH- It was found in Ref. [l(j by an analysis 
of the temperature-dependent Shubnikov-de Haas oscil- 
lations that the strongly increased mass is also the case 
in a dilute 2D electron system in (111) silicon with rel- 
atively high disorder, as indicated by both the absence 
of metallic temperature dependence of zero-field resis- 
tance and the relatively low mobility. Moreover, it was 
found that the relative mass enhancement as a function 
of r s is system- and disorder-independent being deter- 
mined by electron-electron interactions only. However, 
no sign of the enhanced spin susceptibility, which is 
proportional to the product of the g factor and effective 
mass, m, was observed in parallel-field magnetoresistance 
measurements [7J. Rather than tending to zero at a fi- 
nite electron density like in (lOO)-silicon MOSFETs U, 
the resistance saturation field, i? sa t, corresponding to the 



onset of full polarization of the electrons' spins was found 
to increase approximately proportionally with the elec- 
tron density. It is worth noting that the valley degen- 
eracy in (lll)-silicon MOSFETs is equal to g v = 2 (as 
well as the spin degeneracy g s — 2), rather than g v = 6 
(see, e.g., Ref. [llj). as inferred from the filling factor 
change for the period of Shubnikov-de Haas oscillations 
in weak magnetic fields, equal to 4. The reduced valley 
degeneracy was attributed, in particular, to the existence 
of inhomogeneous strains at the interface [I2J. The up- 
per electron subbands can be important for the behav- 
ior of B sat . Filling of these subbands has been observed 
recently in studies of Shubnikov-de Haas oscillations in 
(lll)-silicon samples of a special design: a 2D electron 
system forms at the interface between (111) silicon and 
vacuum [9(. 

In this paper, we study low-temperature parallel-field 
magnetotransport in a 2D electron system in (111) sili- 
con in a wide range of electron densities. We find that 
the polarization magnetic field, B x , obtained by scal- 
ing the weak-parallel-field magnetoresistance at different 
electron densities, tends to zero at a finite electron den- 
sity in a linear fashion. It corresponds to the spin sus- 
ceptibility x °c 9 m that increases strongly at low densi- 
ties, B x oc n s /Xi which is consistent with the increase of 
to observed earlier in this electron system. The polar- 
ization field i? sa t, determined by resistance saturation, 
turns out to deviate to lower values than B x with in- 
creasing electron density. This can be explained by the 
parallel-field-induced filling of the upper electron sub- 
bands corresponding to another two valleys in (111) sili- 
con. The subbands' splitting in our samples is estimated 
at A w 2 meV. 

Measurements were made in an Oxford dilution re- 
frigerator with a base temperature of w 30 mK 
(Chernogolovka) and a dilution refrigerator with a base 
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temperature of » 80 mK (Toulouse) in magnetic fields 
up to 14 T on (lll)-silicon MOSFETs similar to those 
previously used in Ref. Q • Samples had the Hall bar ge- 
ometry with width 400 /im equal to the distance between 
the potential probes. Application of a dc voltage to the 
gate relative to the contacts allowed one to control the 
electron density. Oxide thickness was equal to 154 nm. 
In highest-mobility samples, the normal of the sample 
surface was tilted from [111]- toward [110]-direction by a 
small angle of 8°. Anisotropy for electron transport in 
such samples does not exceed 5% at n s = 3 x 10 11 cm~ 2 
and increases weakly with electron density, staying be- 
low 25% at n s = 3 x 10 12 cm~ 2 [i"o| . The resistance was 
measured by a standard 4-terminal technique at a low 
frequency (5-10 Hz) to minimize the out-of-phase signal. 
Excitation current was kept low enough (below 2 nA) 
to ensure that measurements were taken in the linear 
regime of response. The current was parallel to the in- 
plane magnetic field. Electron densities were determined 
by Shubnikov-de Haas oscillations in perpendicular mag- 
netic fields, for which the sample orientation was changed 
at room temperature and another run was made. The 
threshold voltage was checked to be the same in different 
runs, although the peak electron mobility at T m 1.5 K 
varied in the range between 2500 and 4000 cm 2 /Vs. Two 
samples were used in the experiments; below, we show 
data obtained on one of them. The results reported are 
independent of a sample and a sample state and are re- 
producible in different cryostats. 

Additional measurements were made in a rotator- 
equipped He 4 cryostat at a temperature of « 1.5 K in 
pulsed magnetic fields (Toulouse). In each pulse, the 
magnetic field swept up to « 48 T (or lower) with rising 
and falling times of about 50 ms and 0.3 s, respectively. 
Excitation current with frequency 2 kHz did not exceed 
600 nA. 

Typical curves of the low-temperature magnetoresis- 
tance p(B) in a parallel magnetic field are displayed in 
Fig. [lja). The resistivity increases with field and satu- 
rates above a certain density-dependent magnetic field. 
Deep in the metallic regime (p <C h/e 2 ), the magnetore- 
sistance depends weakly on temperature in the experi- 
mental range of temperatures (at T < 1.5 K). This limit 
is realized for progressively narrower initial intervals on 
the curve p(B), as the electron density is lowered. The 
data discussed in this paper are obtained deep in the 
metallic regime. 

In Fig. Hlb), we show how the normalized weak- field 
magnetoresistance, measured at different electron den- 
sities, collapses onto a single curve when plotted as a 
function of B/B x . The scaling parameter B x has been 
normalized to correspond to the magnetic field at which 
the magnetoresistance at the lowest electron densities in 
the good metallic regime saturates and full spin polariza- 
tion of the electrons is reached [H, [3, EH • Within the 
accuracy with which the saturation field B sat can be de- 
termined, we use for the normalization the curve p(B) at 
h s k6x 10 11 cm~ 2 , measured in pulsed magnetic fields; 
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FIG. 1: (a) Magnetoresistance in dc parallel magnetic fields 
at a temperature of « 80 mK at different electron densities in 
the metallic regime, (b) Scaled curves of the normalized mag- 
netoresistance at different n s versus B / B x for dc (upper set) 
and pulsed (lower set) parallel magnetic fields. The electron 
densities are indicated in units of 10 11 cm -2 . 
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FIG. 2: The normalized magnetoresistance as a function of 
B/Bsat at different electron densities for pulsed parallel mag- 
netic fields. Also shown is the way to determine the saturation 
field. 



n s (10 11 cm" 2 ) 

FIG. 3: Dependence of the fields _B sat (open circles) and B x 
(squares, triangles, and dots) on electron density. Squares and 
triangles are obtained in two runs in dc magnetic fields, and 
dots and open circles correspond to pulsed magnetic fields. 
The dotted lines are linear fits. The fit for B x extrapolates 
to a critical density n x . Also shown by a dashed line is the 
data for the enhanced effective mass of Ref. [10( recalculated 
into B x using the value of g — 2 in bulk silicon. 



we have verified that the fields B allt and B x are practi- 
cally temperature-independent in the range of tempera- 
tures used. At high electron densities, the value of B sat 
is noticeably smaller than B x in this electron system, as 
is evident from Fig. [2] which shows the normalized mag- 
netoresistance versus B/B sat in pulsed magnetic fields. 
Apparently, the scaled data in Fig. [Hb) are a function 
of the degree of spin polarization, £ = B/B x , in weak 
magnetic fields but not at £ » 1. 

In Fig. [3l we show the so-determined polarization fields 
B sa t and B x as a function of electron density. The field 
B x obtained by scaling the weak-field magnetoresistance 
depends linearly on n s and tends to vanish at a finite 
electron density, n x . Our procedure provides good ac- 
curacy for determining the behavior of B x with electron 
density, i.e., the functional form of B x (n s ), even though 
the absolute value of B x is determined not so accurately. 
Based on the form of the spin polarization parameter 
£ = X-B/MB n s = 9I j -bB/2Ef = B/B x (where fig is the 
Bohr magneton and Ep = Trh 2 n s /2m is the Fermi en- 
ergy in B = 0), we have recalculated the data for the 
enhanced effective mass m(n s ) of Ref. [l(| into B x (n s ) 
(dashed line), using the value of g = 2 in bulk silicon 
which is close to the g factor measured in silicon MOS- 
FETs [4]. Both dependences B x (n s ) are in agreement 
with each other, and the spin susceptibility \ <x gm is 
strongly enhanced at low densities. The resistance sat- 
uration field -B sa t deviates to lower values than B x with 
increasing electron density, B sa ±(n s ) being a linear de- 
pendence with approximately twice as small a slope. We 
note that our data for B saii are concurrent with earlier 
results Q • 

The observed behavior of the polarization field B x with 



FIG. 4: The schematic behavior of the fields B sat and B x with 
electron density, as discussed in text. There are break points 
in the dependences B S at(n a ) at Ef = A/2 (or Ef (Bsat) = A) 
and B x (n s ) at Ef = A. 



electron density, which corresponds to the sharply in- 
creased spin susceptibility at low densities, is very simi- 
lar to that in a 2D electron system in (lOO)-silicon MOS- 
FETs [HI]. The same interaction strength in (lll)-silicon 
MOSFETs is expected to be achieved at electron densi- 
ties about four times higher than in (100) silicon, since 
in the first case the effective mass m& = 0.358m e (where 
m e is the free electron mass) in bulk silicon is approx- 
imately twice as large. Based on the critical density, 
n x ~ 8 x 10 10 cm~ 2 , of vanishing polarization field for 
(lOO)-silicon MOSFETs, one expects the critical region 
in (111) silicon at electron densities about 3 x 10 11 cm~ 2 . 
This is consistent with our results for B x (Fig. [3]). 

It is interesting that the polarization field B aa t{n s ) be- 
haves differently from B x (n s ). This finding is in con- 
trast to the case of (lOO)-silicon MOSFETs where the 
fields B x and B sat are practically the same [|| and the 
scaled data are described reasonably well by the theoret- 
ical dependence of the normalized magnetoresistance on 
the degree of spin polarization [l6| in a wide range of £. 
Clearly, the discrepancy between both polarization fields 
in (lll)-silicon MOSFETs reflects the fact that the spin 
polarization parameter is not linear in B at large £ « 1. 
Below, we discuss possible mechanisms for the effect. 

A simple reason for the value of saturation field to 
become reduced as the electron density is increased might 
be the orbital effects that occur when the thickness of the 
2D electron system becomes comparable to the magnetic 
length, l B = {Hc/eB) 1 / 2 . These are expected to lead to 
an increase of the effective mass with parallel magnetic 
field (see, e.g., Ref. [l7|). However, no such reduction in 
the .B sa t(n s ) dependence has been found for (lOO)-silicon 
MOSFETs in the same range of pulsed parallel magnetic 
fields [HI]. Because the corresponding thicknesses of the 
2D electron systems [l9| for both orientations in silicon 
MOSFETs are close to each other, the orbital effects are 
likely to be small for our case. 

Another reason for the observed reduction of B sat is 
the filling of the upper electron subbands in the fully 
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spin-polarized regime. The Fermi energy should double 
at the onset of complete spin polarization [lq |. if Ep is 
below the mid-gap, A/2 (here A is the subband sepa- 
ration). In the opposite case Ep > A/2 the parallel 
magnetic field promotes filling of the upper subbands. 
These are empty in the experimental range of n s in our 
samples in B = (i.e., Ep < A) but can be filled at 
high electron densities in parallel magnetic fields pro- 
vided the condition Ep > A/ 2 is met. In this case the 
Fermi energy at the onset of complete spin polarization, 
Ep{B sat ) = .g/i_B-Bsat, exceeds the splitting A and should 
increase with electron density with a slope reduced by 
the factor of two (see Fig. [4J, as can be expected from 
the results of Ref. 0. There, the filling factor change 
for the period of the weak-field Shubnikov-de Haas os- 
cillations in (lll)-silicon samples has been found to be 
equal to 8, which immediately indicates that the upper 
subbands in question correspond to another two valleys. 
The slope reduction expected for both Ep(B sat ) and B sat 
is in agreement with the data in Fig. [3] Note that the 
observed linear increase of B aat with n s indicates that 
a possible influence of the diamagnetic contribution to 
A is small. Thus, the parallel-field-induced filling of the 
upper subbands corresponding to another two valleys in 
(111) silicon can account for the experimental behavior 
of the saturation field. 

We estimate the splitting of the subbands at A « 
2 meV for n s ps 6 x 10 11 cm~ 2 . This value is ap- 
preciably bigger than the splitting A w 0.6 meV ob- 
tained in Ref. [9( , which is consistent with the pattern of 
Shubnikov-de Haas oscillations in the two sets of samples. 
The fact that the splitting A is not universal in different 
samples may be due to different inhomogeneous strains at 



the interface [12j | which are sensitive to a particular pro- 
cedure of sample fabrication. Based on the experimental 
data and ignoring the possible (linear) dependence of A 
on n s , one can expect filling of the upper subbands in 
B = at electron densities above ks 1.5 x 10 12 cm -2 . 
Once no double the periodicity of Shubnikov-de Haas os- 
cillations in weak magnetic fields was observed in our 
samples up to n s w 2.5 x 10 12 cm -2 yet higher elec- 
tron densities should be tried to clarify the point. Note 
that at very high electron densities about 1 x 10 13 cm -2 , 
observation of the sixfold valley degeneracy was reported 
[20| , although the high g v can be promoted by additional 
annealing of the samples (2~j| . 

In summary, we have found that the polarization field 
B x , obtained by scaling the weak-parallel-field magne- 
toresistance at different electron densities in (lll)-silicon 
MOSFETs, tends to zero at a finite electron density in 
a linear fashion. It corresponds to the spin susceptibility 
that increases strongly at low densities. This is consistent 
with the increase of the effective mass observed earlier 
in this electron system. The polarization field -B sa t, de- 
termined by resistance saturation, turns out to deviate 
to lower values than B x with increasing electron den- 
sity. The reduction of B sa ± can be explained by filling of 
the upper electron subbands in the fully spin-polarized 
regime. The results obtained allow us to estimate the 
subbands' splitting in our samples. 
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